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ABSTRACT 

We examine the constraints on the epoch of reionization (redshift Zj.) set by 
recent WMAP-3 observations of Te = 0.09 ± 0.03, the electron-scattering optical 
depth of the cosmic microwave background (CMB), combined with models of high- 
redshift galaxy and black hole formation. Standard interpretation begins with 
the computed optical depth, Tg = 0.042 ± 0.003, for a fully ionized medium out to 
Zr = 6.1±0.15, including ionized helium, which recombines at z ~ 3. At z > Zr, one 
must also consider scattering off electrons produced by from early black holes (X- 
ray pre-ionization) and from residual electrons left from incomplete recombination. 
Inaccuracies in computing the ionization history, Xe{z) add systematic sources of 
uncertainty in Te. The required scattering at z > Zr can be used to constrain the 
ionizing contributions of "first light" sources. In high- 2; galaxies, the star-formation 
efficiency, the rate of ionizing photon production, and the photon escape fraction are 
limited to producing no more than At^ < 0.03 ±0.03. The contribution of minihalo 
star formation and black-hole X-ray preionization at 2; = 10-20 are suppressed by 
factors of 5-10 compared to recent models. Both the CMB optical depth and H I 
(Lya) absorption in quasar spectra are consistent with an H I reionization epoch 
at Zr ^ Q providing ~ 50% of the total at z < Zr, preceded by a partially ionized 
medium at 2; ~ 6-20. 

Subject headings: cosmology: theory — cosmic microwave background — inter- 
galactic medium 
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1. INTRODUCTION 

In recent years, an enormous amount of exciting cosmological data have appeared, ac- 
companied by theoretical statements about early galaxy formation and the first massive stars. 
Many of these statements were reactions to first-year (WMAP-1) results (Kogut et al. 2003; 
Spergel et al. 2003) from the Wilkinson Microwave Anisotropy Probe (WMAP). These papers 
inferred a high optical depth to the cosmic microwave background (CMB) and suggested early 
reionization of the intergalactic medium (IGM). Other conclusions came from simplified mod- 
els and assumptions about the stellar initial mass function (IMF), atomic/molecular physics, 
radiative processes, and prescriptions for star formation rates and escape of photoionizing 
radiation from protogalaxies. 

In this paper, we focus on the reionization epoch, defined as the redshift Zr when the IGM 

becomes nearly fully ionized over most of its volume (Gnedin 2000, 2004). Our knowledge 
about reionization comes primarily from two types of observations: H I (Lya) absorption in 
the IGM and optical depth of the CMB. Spectroscopic studies of the "Gunn-Peterson" (Lya) 
absorption toward high-redshift quasars and galaxies imply that H I reionization occurred not 
far beyond z ^ 6 (Becker et al. 2001; Fan et al. 2002, 2006), and that He II reionization 
occurred at z ~ 3 (Kriss et al. 2001; ShuU et al. 2004; Zheng et al. 2004). Third-year data 
from WMAP (Spergel et al. 2006; Page et al. 2006) suggest that reionization might occur at 
z 10, although large uncertainties remain in modehng of the CMB optical depth. 

The WMAP and Lya absorption results are not necessarily inconsistent, since they probe 
small amounts of ionized and neutral gas, respectively. In addition, both the H I absorbers and 
ionized filaments in the "cosmic web" (Cen & Ostriker 1999) are highly structured at redshifts 
z < 10 and affect the optical depths in Lya. For example, in order to effectively absorb all 
the Lya radiation at ^ ~ 6 requires a volume- averaged neutral fraction of just xhi ~ 4 x 10~^ 
(Fan et al. 2006). Simulations of the reionization process (Gnedin 2004; Gnedin & Fan 2006) 
show that the transition from neutral to ionized is extended in time between 2; = 5 — 10. The 
first stage (pre-overlap) involves the development and expansion of the first isolated ionizing 
sources. The second stage marks the overlap of the ionization fronts and the disappearance of 
the last vestiges of low-density neutral gas. Finally, in the post-overlap stage, the remaining 
high-density gas is photoionized. The final drop in neutral fraction and increase in photon 
mean free path occur quite rapidly, over Az ~ 0.3. 

Initial interpretations of the high WMAP-1 optical depth were based on models of "sudden 
reionization", neglecting scattering from a partially ionized IGM ect z > Zr- The revised optical 
depths, Te = (0.09 - 0.10) ± 0.03 (Spergel et al. 2006; Page et al. 2006) also assume a single 
step to complete ionization (xg = 1 at redshift Zr), although Spergel et al. (2006) explored a 
two-step ionization history, resulting in a broader distribution of Tg. However, the curves 
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for two-step ionization provide little constraint on the redshift of reionization, or on an IGM 
with low partial ionization fractions, <S 0.1, as we discuss in § 3. Figure 3 of Spergel et al. 
(2006) shows that the parameters Xe and Zr are somewhat degenerate, each with a long tail 
in the likelihood curves. This emphasizes the possible importance of contributions to Tg from 
ionizing UV photons dX z > Zr from early massive stars and X-rays from accreting black holes 
(Ricotti & Ostriker 2004; Begelman, Volonteri, & Rees 2006). 

Models of the extended recombination epoch (Seager, Sasselov, & Scott 2000) predict a 
partially ionized medium at high redshifts, owing to residual electrons left from incomplete 
recombination. These electrons produce additional scattering, Arg ~ 0.06 {z 10-700). 
These effects must be computed in CMB radiation transfer codes such as CMBFAST and 
RECFAST, but only a portion of this scattering affects the large angular scales {I < 10) where 
WMAP detects a polarization signal. X-ray preionization can also produce CMB optical depths 
Te > 0.01 (Venkatesan, Giroux, & ShuU 2001, hereafter VGSOl; Ricotti, Ostriker, & Gnedin 



In § 2.1, we show that a fully ionized IGM from 2; = to the reionization epoch at 
Zr = 6.1 ± 0.15 (Gnedin & Fan 2006), produces optical depth, Tg ^ 0.042 ± 0.003, nearly half 
the WMAP-3 value. Therefore, the high-redshift ionizing sources are hmited to producing an 
additional optical depth, Arg < 0.03 ±0.03. In § 3, we discuss the resulting constraints on the 
amount of star formation and X-ray activity at 2; > 7 and limits on star formation in minihalos. 
The IGM probably has a complex reionization history, with periods of extended reionization 
for H I and He II (Venkatesan, Tumlinson, & Shull 2003; Gen 2003; Wyithc & Locb 2003; Hui 
& Haiman 2003). Because Tg measures the integrated column density of electrons, there are 
many possible scenarios consistent with the current (WMAP-3) level of GMB data. 



To elucidate the dependence of GMB optical depth on the epoch of reionization (redshift 
Zr)i we integrate the electron scattering optical depth, Te{zr)i for a homogeneous, fully ionized 
medium out to Zr- For instantaneous, complete ionization at redshift z^, we calculate Tg as the 
integral of necrrdi, the electron density times the Thomson cross section along proper length. 



^0 

We adopt a standard AGDM cosmology, in which {d£/dz) — c{dt/dz) = (1 -|- z)''^[c/ H{z)], 
where H{z) — ifo[^^m(l + + ^a]^''^ and + Q\ — 1 (no curvature). The densities of 



2005). 



2. OPTICAL DEPTH TO ELECTRON SCATTERING 



2.1. Analytic Calculation of Optical Depth 




(1) 
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hydrogen, helium, and electrons are written uh — [ii—Y)Pcr/'>TiH] i^+z)^, n^B.e — ynn, and Ug — 
nnil + y), if helium is singly ionized. We also add the small {t^ ~ 0.0020) contribution from 
doubly- ionized helium (He 111) at rcdshifts z < 3 (Shull et al. 2004). We assume a helium mass 
fraction by mass Y = 0.244 and define y = {Y/A)/{1 — Y)^ 0.0807 (He fraction by number). 
The critical density is pcr = (1-879 x 10"^^ g cm-^)h'^ where h = (i^o/lOO km s"^ Mpc"^). The 
above integral can be done analytically: 

= (i) (^) [ """Z'^'l . P) 

where {1 — Y){1 + y) = (1 — 3F/4). For large redshifts, 0^,(1 + z)^ 3> ^a, and the integral 
simplifies to 
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Here, we substituted the WMAP-3 parameters (Spergel et al. 2006): f^b/i^ ^ 0.0223 ± 0.0009, 
^mh^ — 0.1271q ;o?3, and h — 0.73 ± 0.03, and adopted a primordial helium abundance Yp — 
0.244 ± 0.002 (Izotov & Thuan 1998; Olive & Skillman 2001). A recent paper (Peimbert et al. 
2007) suggests a revised Yp — 0.2474 ± 0.0028 based on new He I atomic data. 

Prom the approximate expression (eq. 3), we see that Tg oc (pcr^b^m Hq ). Thus, Tg 
is nearly independent of the Hubble constant, since Pcr oc h? while the combined parameters, 
Q^bh^ and O^mh^, are inferred from D/H, CMB, and galaxy dynamics. The scahng with h 
cancels to lowest order; a slight dependence remains from the small Qa term in equation (2). 

If we invert the approximate equation (3), we can estimate the primary reionization red- 
shift, (1 -\-Zr) ~ (6.6)[re(2;r)/0.04]^/^, where we scaled to the value, Tg = 0.04, expected for full 
ionization back to 2;^ ~ 6. As discussed in § 2.2, this is approximately the WMAP-3 value of 
optical depth {je ~ 0.09) reduced by Arg ~ 0.05. This extra scattering may arise from high-z 
star formation. X-ray preionization, and residual electrons left after incomplete recombination. 
The latter electrons are computed to have fractional ionization x^. (0.5 — 3.0) x 10~^ be- 
tween z = 10-700 (Seager et al. 2000). Inaccuracies in computing their contribution therefore 
add systematic uncertainty to the CMB-derived value of Tg. Partial ionization may also arise 
from the first stars (Venkatesan, Tumlinson, & Shull 2003, hereafter VTS03) and from pene- 
trating X-rays produced by early black holes (VGSOl; Ricotti & Ostriker 2004, 2005). These 
additional ionization sources contribute electron scattering that must be subtracted from the 
WMAP-3 values. They will lower the reionization redshift, Zr — IO.7II3, derived (Spergel et 
al. 2006) in the absence of partial ionization at z > Zr, and they may bring the WMAP-3 and 
Gunn-Peterson results into agreement for the epoch of complete reionization. 

In our calculations, described in § 3, we make several key assumptions. First, we assume 
a fully ionized IGM out to Zr ~ 6, accounting for both H+ and ionized hehum. (Helium 
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contributes 8% to Te, assuming He II at 2; > 3 and an additional Tg ~ 0.002 for He HI at 
2; < 3). Second, we investigate the effects of IGM partial ionization at 2; > Zr- Finally, in 
computing the contribution of residual electrons at high redshifts, we adopt the concordance 
parameters from the WMAP-3 data set. The CMB optical depth is formally only a Scr result, 
which may change, as WMAP produces better determinations of the matter density, Qm, 
and the parameters, as and n^, that govern small-scale power. Both as and n have well- 
known degeneracies with Tg in CMB parameter extraction (Spergel et al. 2006). Therefore, 
their derived values may change in future CMB data analyses, especially as the constraints 
on Te continue to evolve. In addition, inaccuracies in the incomplete recombination epoch 
and residual ionization history, Xe{z), add uncertainties to the CMB radiative transfer, the 
damping of £-modes, and the polarization signal used to derive an overall Tg. 



2.2. Residual Electrons in the IGM 

We now discuss the contribution of residual electrons in the IGM following the recombina- 
tion epoch at z ^ 1000. Scattering from these electrons is significant and is normally accounted 
for in CMB transport codes such as CMBFAST (Seljak & Zaldarriaga 1996) through the post- 
recombination IGM ionization history, Xe{z). However, a number of past papers are vague on 
how the ionization history is treated, which has led to confusion in how much residual optical 
depth and power-damping has been subtracted from the CMB signal. Modern calculations 
of how the IGM became neutral have been done by Seager et al. (2000), although their code 
(RECFAST) continues to be modified to deal with subtle effects of the recombination epoch 
and the atomic physics of hydrogen (2s — > Is) two-photon transitions (W. Wong & D. Scott, 
private communication) . 

To illustrate the potential effects of high-2 residual electrons, we have used numbers 
from Figure 2 of Seager et al. (2000), the top-panel model, which assumed a cosmology with 
fttot = 1, = 0.05, h = 0.5, Y = 0.24, and Tcmb = 2.728 K. At low redshifts, z ^ Zr, just 
before reionization, they find a residual electron fraction Xe,o ~ 10"^-^. We fitted their curve 
for log Xe out to 500 to the formula: 

Xe{z) = Xe,o 10°-°'^^(^+^) ^ {5 X 10"^) exp[a(l + z)] , (4) 

where a ^ 2.303 x 10~^. More recent recombination calculations (W. Y. Wong & D. Scott, 
private communication) using WMAP-3 parameters (Qb = 0.04, h = 0.73, flrn — 0.24, Qa = 
0.76, Y = 0.244) find somewhat lower values, Xefl ~ lO"^-^^ with a ^ 2.12 x 10"^. We 
attribute the lower Xe^ to the faster recombination rates arising from their higher assumed 
baryon density, Qbh'^ = 0.0213, compared to Qbh'^ = 0.0125 in Seager et al. (2000). 
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To compute the electron-scattering of the CMB from these "frozen-out" electrons, we use 
the same integrated optical depth formula (equation 1), in the high-^ limit, where {M/dz) ~ 
i (c/i/o)^m ^^(1 + We integrate over the residual-electron history. 



[\^z)-\clR{z)\ 

7 back to a final redshift Zj ^ z^.^ io find 



from Zr 



(Are), 
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1 zfl"^ Xe{z) dz 



ril+zf) 

(3.27 X 10"^) Xe,o / u^'^ exp{au) du . 



(5) 



A rough estimate to the residual scattering comes from setting a ~ and adopting a constant 
ionized fraction Xe{z), 



(ATe)res ~ (2.36 X 10"^) [(1 + Zff/^ - (1 + Zrf/^]{Xe 



(6) 



This estimate gives Tg = 0.044 for Zr — 6, Zf = 700, and (xg) ~ 10~^. More precise values of 
Tg can be derived from the exact integral (eq. 5) by expanding the exponential as a sum and 
adopting the limit Zf ^ Zr, 



(ATe)res = (0.0184) 



n + 3/2) n! 



501 



E 

n=0 



(7) 



Prom Zr — 7 out to various final redshifts, we find cumulative optical depths: (ATe)res = 
0.0041 (from 7 < z < 200), (At^),^, = 0.0088 (from 7 < z < 300), (Are)res = 0.0157 (from 
7 < z < 400), and (ATe)res = 0.0256 (from 7 < z < 500). For z > 500, the approximate 
formula (eq. 4) underestimates Xe, but one can integrate the appropriate curves (Seager et al. 
2000) using piecewise-continuous linear fits. Between 500 < z < 600, wc find Xe = (2.94 x 
10-^)exp[0.003224(l + z)], and for 600 < ;2 < 700, x^ = (1.28 x 10"^) exp[0.004606(l + z)]. 
Integration then yields additional contributions of Arg ^ 0.0135 for z — 500-600 and Are ~ 
0.021 for z — 600-700. These calculations therefore give a total optical depth in residual 
electrons Te ~ 0.06 back to 2; = 700. These electrons have maximum influence on angular 
scales with harmonic ^max ~ 2^;^/^ 20-50 (Zaldarriaga 1997). At higher redshifts, Xe rises to 
10-2-^ at ^ = 800 and to lO"^-^ at ^ = 1000, where the CMB source function will affect the 
"free-streaming" assumption used in CMBFAST (Seljak & Zaldarriaga 1996). 



3. IMPLICATIONS FOR REIONIZATION MODELS 



The WMAP-3 measurements of fluctuations in temperature (T) and polarization (E) have 
been interpreted to estimate totoZ electron-scattering optical depths of Tg = (0.09 — 0.10) ±0.03. 
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The central values, Te = 0.09 (Spergel et al. 2006) and Tg = 0.10 (Page et al. 2006) come, 
respectively, from computing the likelihood function for the six-parameter fit to all WMAP 
data (TT, TE, EE) and for just the EE data as a function of Te- Because of the challenges in 
translating a single parameter (rg) into a reionization history, Xe{z), it is important to recognize 
the sizable error bars on Tg. At the 68% confidence level, Tg could range from "low values" 
(0.06-0.07) up to values as high as 0.12-0.13. 

In § 2.1, we showed that ~ 50% of this Tg can be accounted for by a fully ionized IGM 
at z < Zr- Recent Gunn-Peterson observations of 19 quasars between 5.7 < z < 6.4 (Fan 
et al. 2006; Gncdin & Fan 2006) arc consistent with a reionization epoch of 2;^ = 6.1 ± 0.15. 
According to equation (2), this produces Tg = 0.042 ± 0.003, where our error propagation 
includes relative uncertainties in Zr (2.5%), (8.4%), and fife/i^ (4.0%). Residual post- 

recombination electrons produce a substantial optical depth from z 10 back to z ^ 700, 
which uniformly damps all angular scales. However, their effect on the TE and EE power is 
considerably less on large angular scales (£ < 10). Thus, we can characterize a portion of 
the WMAP-3 observed optical depth, Tg ~ 0.09 ± 0.03, through known sources of ionization. 
The "visible ionized universe" out to Zr = 6.1 accounts for Tg = 0.042, while high-z partial 
ionization could contribute anywhere from Tg ^ 0.01 to 0.06. We therefore assume that an 
additional optical depth, Arg < 0.03 ± 0.03, can be attributed to star formation and early 
black hole accretion at z > Zr- 

Our calculations represent an important change in the derivation of z.,. from Tg, suggesting 
that the amount and efficiency of high- 2; star formation need to be suppressed. This suggestion 
is ironic, since WMAP-1 data initially found a high Tg = 0.17 ± 0.04 (Spergel et al. 2003) 
implying a surprisingly large redshift for early reionization, ranging from 11 < 2;^ < 30 at 
95% confidence (Kogut et al. 2003). These results precipitated many investigations of star 
formation at 2; = 10 — 30, some of which invoked anomalous mass functions, very massive 
stars (VMS, with M > 140 Mq), and an increased ionizing efficiency from zero-metallicity 
stars (VTS03; Wyithe & Loeb 2003; Gen 2003; Giardi, Ferrara, & White 2003; Sokasian et al. 
2003, 2004). Tumlinson, Venkatesan, & Shull (TVS04) disputed the hypothesis that the first 
stars had to be VMS. They showed that an IMF dominated by 10 - 100 Mq stars can produce 
the same ionizing photon budget as VMS, generate GMB optical depths of 9-14%, and still 
be consistent with nucleosynthetic evidence from extremely metal-poor halo stars (Umeda & 
Nomoto 2003; Tumlinson 2006; Venkatesan 2006). 

Although the IGM recombination history, Xe{z), is included in calculations of GMBFAST 
and in GMB parameter estimation, the best-fit values of Tg from WMAP-3 and earlier GMB 
experiments have been attributed exclusively to the contribution from the first stars and/or 
black holes at 2; < 20. The contributions from post- recombination electrons (20 < z < 1100) 
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have not always been subtracted from the data. This post-recombination contribution was 
relatively small in some earlier models of reionization (Zaldarriaga 1997; Tegmark & Silk 
1995) that explored optical depths of Te = 0.5 — 1.0 and suggested reionization epochs up to 
~ 100. However, with current data indicating late reionization, it becomes particularly 
important to consider contributions to Tg prior to the first sources of light. 

The new WMAP-3 results find a lower Tg, but they also suggest less small-scale power 
available for reionizing sources, owing to lower normalization parameters, cig 0.74lQQg and 
flmh^ ~ 0.1271;q q53- This reduction is somewhat offset by the reduction in spectrum tilt from 
ris = 0.99 ± 0.04 (WMAP-1) to = 0.95ll[^:J]J^ (WMAP-3). Alvarez et al. (2006) argue, 
from the lower values of and erg, that both WMAP-3 and WMAP-1 data require similar 
(high) stellar ionizing efficiencies. Haiman & Holder (2003) use the lower Tg to suggest that 
massive star formation was suppressed in minihalos. Our results on a lower Are make these 
requirements even more stringent, as we now quantify. 

Semi-analytic and numerical models of reionization (Ricotti, Gnedin, & ShuU 2002a,b; 
VTS03, Haiman & Holder 2003) show that the efficiency of ionizing photon injection into the 
IGM can be parameterized by the "triple product", A^-y/*/esc- Here, /* represents the star- 
formation efficiency (the fraction of a halo's baryons that go into stars), Nj is the number 
of ionizing photons produced per baryon of star formation, and /esc is the fraction of these 
ionizing photons that escape from the halo into the IGM. We can now use our calculations to 
constrain the amoimt of high- 2; star formation through the product of these three parameters, 
henceforth referred to as the "efficiency". For the ionization history in equation (5), we set 
Xe{z) = N^f^fcscCL{,z)fh{z), where cl{z) is the space-averaged baryon clumping factor of ionized 
hydrogen, cl = (^mi) (^hii)^- We assume that cl is the same for H H and He HI. The factor 
fi,{z) is the fraction of baryons in collapsed halos, computed through the Press-Schechter 
formalism (as in VTS03) for the cosmological parameters from WMAP-3. We assume that 
Njf^fesc is constant with redshift. 

There is surely some dependence of each of these parameters on the halo mass and envi- 
ronment (Haiman & Bryan 2006; Ricotti & ShuU 2000). Since we have already parameterized 
the intrahalo recombinations through /esc, we account for the loss of ionizing photons on 

IGM scales through Cx, in two forms: (1) a power-law form with slope (3 = —2 from the semi- 
analytic work of Haiman & Bryan (2006); and (2) the numerical simulations of Kohler, Gnedin 
& Hamilton (2006), using their case C (overdensity 5 ~ 1 for the large-scale IGM) for Cr, the 
recombination clumping factor corresponding to our definition of C/,. In the latter case, the 
clumping factor is almost constant {cl ~ 6) until the very end of reionization. Together, these 
two different cases provide bounds on the range of possible values. 

With these assumptions, we can use equation (5) and the allowed additional optical depth. 



- 9 - 



Are < 0.03 ± 0.03, to constrain the ionizing efficiency of tfie first stars. In Figure 1, we plot 
tlie efficiency as a function of At^, for star formation in Lya-cooled lialos (virial temperature 
Tvir > IC^ K) and in H2-cooled minilialos (T^ir > 10^ K). We consider the two clumping 
factors from Kohler et al. (2006) and Haiman & Bryan (2006). For j\ = f^sc = 0.1, we 
indicate the efficiencies corresponding to two cases of interest, assuming a Salpeter initial mass 
function (IMF): (1) A^^ = 60,000 for a metal-free IMF (10-140 Mq) that agrees with both 
CMB and nucleosynthetic data; and (2) A^^^ = 4000 for a present-day IMF (1-100 Mq). Note 
that — 34,000 is consistent with the WMAP-3 inferrence, Te ~ 0.1, if the photons arise 
from zero- metal stars (Tumhnson 2006). These values of Nj were derived (TVS04) from the 
lifetime-integrated ionizing photon production from various stellar populations and IMFs and 
used as inputs in cosmological reionization models. 

Figure 1 shows that fairly modest efficiencies of massive star formation are consistent with 
limiting the ionizing contribution of minihalos to Are < 0.03 — 0.06. This is consistent with 
suppression of star formation and ionizing photon production in mini-halos (Haiman & Bryan 
2006). Much larger efficiencies, close to those of metal-free stellar populations, are required 
for larger halos. Interestingly, it seems to make little difference what form is assumed for the 
clumping factor, Cl, in constraining the ionizing efficiency of the first stars. This is largely a 
statistical effect arising from the insensitivity of the average optical depth, or electron column 
density, averaged over many beams passing through a clumpy medium. 

The CMB optical depth can also constrain the level of X-ray preionization from high- 
redshift black holes. Ricotti et al. (2005) were able to produce large optical depths, Tg ~ 0.17, 
using accreting high- 2; black holes with substantial soft X-ray fluxes. Their three simulations 
(labeled M-PIS, M-SNl, M-SN2) produced hydrogen preionization fractions Xe = 0.1 - 0.6 
between z — 15 and z — 10, with large co-moving rates of star formation, (0.001 — 0.1) Mq 
Mpc"^ yr~^, and baryon fractions, a;BH ~ 10~^ - 10~^ accreted onto black holes. Any significant 
contribution to Tg from X-ray preionization requires Xe > 0.1. Therefore, the lower value 
of Tg from WMAP-3 reduces the allowed X-ray preionization and black-hole accretion rates 
significantly compared to these models. Ricotti et al. (2005) demonstrated that high optical 
depths (Te ~ 0.17) could be achieved by black-hole X-rays. In their models, the IGM at z > Zr 
was highly ionized {xe ^ 0.01). In this limit, there were few X-ray secondary electrons and 
most of the X-ray energy went into heating the ionized medium. By contrast, the WMAP-3 
optical depth suggests that the IGM was much less ionized at 2; > 2;^. 

Relating the effects of X-ray ionization from early black holes to a Arg and a related X-ray 
production efficiency factor is less straightforward compared to the star formation case, for the 
following reasons. First, unlike ionization by UV photons. X-ray ionization is non-equilibrium 
in nature and the timescales for X-ray photoionization at any epoch prior to 2; = 6 typically 
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exceeds the Hubble time at those epochs (VGSOl). Therefore, a one-to-one correspondence 
between X-ray production at an epoch, and the average efficiency of halos at that epoch is more 
difficult to establish relative to the UV photon case. In addition. X-ray ionization (whether 
from stars or black holes), at least initially when Xg < 0.1-0.2, will be dominated by secondary 
ionizations from X-ray-ionized helium electrons rather than from direct photoionization. This 
may therefore constrain the physical conditions in the IGM (e.g., the level of He ionization) 
rather than those in the parent halo, when we attempt to translate a Arg into an X-ray 
ionization efficiency. Thus, it may be difficult to make precise inferences about the black hole 
density and accretion history from Tg. 

We define an efficiency parameter, ex, for X-rays analogous to the previous case for massive 
star formation. Here, ex is the product of the average fraction of baryons in black holes (in 
halos at 2; > 7) and the number of X-ray photons produced (per baryon accreted onto such black 
holes). We assume that the clumping factor (cl) and escape fraction (/esc) are roughly unity 
for X-rays, given their high penetrating power relative to UV photons. Thus, we define Xe{z) ~ 
exfb{z) for X-rays, where electrons come from H+, He"*", and He"*"^. We assume that each X-ray 
photon produces ~12 hydrogen ionizations, primarily through secondary ionizations from X- 
ray photoclcctrons (VGSOl). Figure 2 shows the allowed additional optical depth, analogous 
to the constraints of Figure 1, for X-ray efficiency in both Lya-cooled halos and minihalos. A 
comparison with Figure 1 reveals that X-rays are capable of much higher ionization efficiency 
relative to Pop HI or Pop II star formation. 

In summary, we have shown that the revised (WMAP-3) values of CMB optical depth, 
Te = 0.09 ± 0.03, lead to a more constrained picture of early reionization of the IGM. Approx- 
imately half of the observed Te comes from a fully ionized IGM back to 2;^ = 6.1 ± 0.15. The 
remainder probably arises from the first massive stars and from accretion onto early black holes 
at z > Zr- Some of the observed Tg may come from scattering from residual electrons left from 
recombination; inaccuracies in computing this ionization history add systematic uncertainty 
to the CMB inferred signal. We have assumed extra scattering, Arg = 0.03 ± 0.03 at 2; > Zr, 
and used this to constrain the effciciencics for production and escape of ionizing photons, both 
from the ffist massive stars (Figure 1) and early black holes (Figure 2). 

In both cases, the picture is of a partially ionized IGM at redshifts z — 6 — 20. For X-ray 
pre-ionization by early black holes, equation (6) can be used to provide an estimate of the 
effects of partial reionization. Between redshifts Z2 ~ 20 and Zi ^ 6, an IGM with ionized 
fraction Xefl = 0.1 (Ricotti et al. 2005) would produce (Arg) = (0.018)(a;e,o/0.1), which is a 
significant contribution to the observed Tg = 0.09 ± 0.03. Ricotti et al. (2005) suggested a 
large Xg ~ 0.1 — 0.6 and pushed their black-hole space densities and accretion rates to large 
values in order to reach the WMAP-1 estimates of Tg = 0.17. Because such large values of Tg 
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are no longer required, the black hole densities and IGM ionization fractions are hkely to be 
considerably less. 

All these constraints depend heavily on uncertain parameterizations of the efficiency of 
star formation and ionizing photon production. However, with more precise measurements 
of CMB optical depth from future missions, there is hope that more stringent constraints on 
high- 2; star formation and black-hole accretion will be possible. 
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Fig. 1. — Efficiency factor, N^f^:fesc, for production and escape of photoionizing radiation 
vs. allowed additional optical depth from ffist stars, Are a.t z > Zr- Efficiency is defined as in 
Haiman & Bryan (2006) for star formation in Lya-cooled halos (Tvir > 10"^ K) and in H2-cooled 
minihalos {T^i^ > 10^ K). Solid and dashed curves correspond to clumping factors from Kohler 
et al. (2006) and Haiman & Bryan (2006) respectively. Horizontal solid lines correspond to 
the efficiencies (TVS04) from a metal-free (10-140 Mq) Pop HI IMF and a standard Salpeter 
(1-100 Mq) Pop II IMF, with = /esc = 0.1. See text for discussion. 



- 15 - 




Fig. 2. — Required X-ray production factor, ex, versus electron-scattering optical depth, Are, 
at redshifts z > Zr. This factor is taken to be the baryon fraction in black holes times the 
photons per accreted baryon (see § 3 for details) and controls the rate of producing escaping 
X-rays from halos with virial temperatures of 10^ K and 10'^ K. The electron ionized fraction 
is related to the total baryon fraction by Xe{z) = exfb{z). 



